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The origin of the highest energy cosmic rays remains an enigma. They offer a window to 
new physics, including tests of physical laws relevant to their propagation and interactions, 
at energies unattainable by terrestrial accelerators. They must be accelerated locally, as 
otherwise background radiations would severely suppress the flux of protons and nuclei, at 
energies above the Greisen-Zatsepin-Kuzmin (GZK) limit (~ 60EeV=6 x lO^^eV). Nearby 
Gamma Ray Bursts (GRBs), Hypernovae, Active Galactic Nuclei (AGNs) and their flares, 
have all been suggested and debated as possible sources. A local sub-population of type Ibc 
supernovae (SNe) with mildly reiativistic outflows have been detected as sub-energetic GRBs 
or X-Ray Flashes (XRFs) and recently as radio afterglows without detected GRB counter- 
parts. We measure the size-magnetic field evolution, baryon loading and energetics, using 
the observed radio spectra of SN 2009bb. We place such engine-driven SNe above the Hillas 
line and establish that they can readily explain the post- GZK UHECRs. 

The highest energy cosmic rays pack such a large amount of energy and have such a low 



flu)Pthat direct detection by satellite-borne instruments is infeasible, allowing them to be detectecP 
only by cosmic ray air showers^. UHECRs beyond the GZK limi^^have been invoked to propose 
tests of known physical laws and symmetries^. Understanding their origin is a crucial step in 
using them as probes of new physics. But, the sources of the highest energy cosmic rays pose an 
intriguing problem, since the magnetic rigidity of these particles are such that the magnetic fields 
in our galaxy are neither strong enough to contain them nor bend them sufficiently^. Yet, among 
the UHECRs which have been detected until now, no concentration have been found towards the 
Milky Way. Hence, it is anticipated that their sources are extragalactic. However, UHECR protons 
with energies above 60 EeV can interact with a significant fraction of CMB photons via the A 
resonance. The cross section of this interaction is such that only those extragalactic cosmic ray 
sources locally (within 200 Mpc of the Earth) can contribute significantly to the flux of UHECRs 
above the so called GZK limiP^. At least 61 sources are required by observations of UHECRs 
until now, to provide cosmic rays with energies above the GZK limit^ Since particles of such high 
energy could not have traveled to the Earth from cosmological distances, unless Lorentz invariance 
breaks down at these energies^, they encourage the search for potential cosmic ray accelerators in 
the local Universe. AGNs have been considered^ as UHECR sources. But, most of them are 
not luminous enoughP, leaving proposed very intense, short duration AGN flares, which are yet 
to be observed, as possible sources^. Classical Gamma Ray Bursts (GRBs) are also considered 
as possible source^^, but most of them occur beyond the GZK horizon^ and cannot contribute 
significantly to the local flux beyond the GZK limit^. Hypernovae have also been suggested as 
sources where particles are boosted to successively higher energies in an ejecta profile extending 
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upto mildly relativistic velocitie^'^. However, they require excessive explosion energy and fail 
to reproduce the flat injection spectrum of UHECRs (see Suppl. Info for discussion of some of the 
proposed sources). 

Soon after their suggestion that SNe come from collapse of a normal star to a neutron star, 
Baade and Zwicky went on to suggest that SNe may be the source of cosmic rays as welf^. Since 
then, SNe and SN remnants have been studied as sources of high energy cosmic rays. However, 
ordinary SNe and their remnants can not produce UHECRs due to two fundamental limitations. 
Firstly, they well lie below the line representing the combination of size and magnetic field required 
to confine and accelerate Iron nuclei with energies of 60 EeV, in the so called Hillas diagramP' 
(Figure 1). The second, even more restrictive, condition obviously not fulfilled by ordinary SNe is 
because ordinary SNe have /3/r ~ 0.05 (P = v/c and F = 1/^/1-/3'^, where v is the speed of 
the blastwave and c is the speed of light in vacuum) which restricts the highest energy cosmic rays 
accelerated in ordinary SNe to well below the GZK limit. 

Until recently, SNe with relativistic ejecta have been spotted exclusively through Long GRBs 
associated with them like GRB 98042^or its twin ORB 031203. The discovery of XRF 06021^1 
associated with SN 2006aj showed that mildly relativistic SNe are hundred times less energetic but 
thousand times more common (in their isotropic equivalent rate, which is relevant for UHECRs 
actually reaching the observer) than classical GRB^. Radio follow up of SNe Ibc have now dis- 
covered the presence of an engine driven outflow from SN 2009bb''^, without a detected GRB. The 
mildly relativistic SNe, detected either using XRFs or radio afterglows, a subset of SNe Ibc are far 
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more abundant at low redshifts required for the UHECR sources, than the classical GRBs. More- 
over, given their mildly relativistic nature, they have the most favorable combination of /3/r ~ 1, 
unlike both non-relativistic SNe and ultra-relativistic classical Long GRBs. 

In order to derive the highest energy upto which these relativistic SNe can accelerate cosmic 
rays, we have to determine the evolution of the size and the magnetic field in the blast-wave. It has 
been demonstrated that a Synchrotron Self Absorption (SSA) model fits the initial radio spectrum 
of SN 2009bb rather well,^' with a low frequency turnover defining the spectral peak shifting to 
lower frequency with time, characteristic of the expansion of the shocked region that powers the 
radio emission. This allows us to measure the evolution of the radii and magnetic fields from 
VLA and GMRT data (see Suppl. Info) at 5 epochs, plotted on the Hillas diagram (Figure 1). This 
clearly demonstrates that SN 2009bb and XRF 060218 can both confine UHECRs and accelerate 
them to highest energies seen experimentally. At the time of the earliest radio observation^ with 
its fortunate combination of /3/r ~ 1, SN 2009bb could have accelerated nuclei of atomic number 
Z to an energy of ~ 6.5 x Z EeV. For example, the source could have accelerated protons. Neon, 
and Iron nuclei to 6.4, 64 and 166 EeV respectively. In this scheme, the highest energy particles are 
likely to be nuclei heavier than protons, consistent with the latest results indicating an increasing 
average rest mass of primary UHECRs with energ)^^. Therefore, our results support the claimed 
preference of heavier UHECRs at the highest energies of the Auger collaboration, although this 
claim is disputed by another experimenP^. 

To estimate whether there are enough relativistic SNe to explain the target objects associated 
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with the ~ 60 detected UHECRs, we require the rate of such transients. SNe Ibc occur at a 
rateP"23of ~ 1.7 x 10"^ Gpc^'^ yr^^. The fraction of Ibc which have relativistic outflows is still a 
somewhat uncertain number, estimatecP^ to be around ~ 0.7%. Hence the rate of SN 2009bb-like 
mildly relativistic SNe is ~ 1.2 x 10~^ Mpc~^ yr^^, which is comparable to the rate of mildly 
relativistic SNe detected as sub-energetic GRBs or XRFs of ~ 2.3 x 10^'' Mpc^^ yr^^. This gives 
us ~ 4 (or 0.5) such objects within a distance of 200 (or 100) Mpc every year. Since SN 2009bb is 
still a unique object, only a systematic radio survey can establish their cosmic rate and statistical 
properties (see Suppl. Info). However, cosmic rays of different energies have different travel delays 
due to deflections by magnetic fields. For a conservative mean delay'^of (r^e/aj/) ~ 10^ yrs we may 
receive cosmic rays from any of 4 (or 0.5) x 10^ possible sources at any point in time. Given the 
situation, in which a direct association between a detected UHECR and its source is unlikel)^^, the 
literature in the subject has focused on the constraintsP'^ placed on plausible sources. We have 
shown in our work that indeed this new class of objects satisfy all these constraints. 

Nuclei are also subject to photo-disintegration by interaction with Lorentz boosted CMB 
photons and can travel upto a distance of ~ 100 Mpc (see Suppl. Info), smaller than but comparable 
to the GZK horizon. So, the local rate of mildly relativistic SNe is high enough to provide enough 
(^ 60) independent sources of cosmic rays with energies above the GZK limit. The value of 
{jdeiay) ^Iso impUcs that it will not be possible to detect UHECRs from a known relativistic SN, 
such as SN 2009bb, within human timescales. However, high energy neutrinos from photo-hadron 
interaction at the acceleration site may be a prime focus of future attempts at detecting these sources 
with neutrino observatories like the IceCube (see Suppl. Info). 
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The required energy injection rate per logarithmic interval in UHECRs^ is Tinj = (0.7 — 
20) X lO'''* erg Mpc^-^ yr^^. Given the volumetric rate of mildly relativistic SNe in the local uni- 
verse, if all the energy injected into UHECRs is provided by local mildly relativistic SNe, then 
each of them has to put in around Esn = (0.3 — 9) x 10^^ ergs of energy, which is comparable 
to the kinetic energy in even a normal SN and can easily be supplied by a coUapsar modeP^. The 
minimum energy in the relativistic outflow of SN 2009bb, required to explain the radio emission 
alone, was founcPtobe E^g ^ 10^^ ergs. Moreover, the mildly relativistic outflow of SN 2009bb 
has been undergoing almost free expansion for ~ 1 year. Our measurements of this expansion 
allows us to show (see Suppl. Info) that this relativistic outflow, without a detected GRB, is signif- 
icantly baryon loaded and the energy carried by the relativistic baryons is EBaryons ^ 3.3 x 10^^ 
ergs. 

If a relativistic outflow carries similar energies in protons, electrons and magnetic fields, the 
radiative cooling of the electrons will lead to an X-ray transienP^. In our model, the acceleration 
occurs in the forward shock produced by the engine driven relativistic ejecta, characterized by a 
single bulk Lorentz factor. Protons and nuclei in such a coUisionless shock show a flat spectrum 
of UHECRp^, consistent with the extragalactic component of the cosmic ray spectrum^"*. Our 
radio observations of SN 2009bb constrain the energy carried by the radiating electron^ and the 
energy of the relativistic baryons powering the almost free expansion for ~ 1 year until now (see 
Suppl. Info). These observations indicate a spectral index of ~ 1 in the optically thin part of the 
radio spectrum^^. This implies a power law distribution of relativistic electrons, with an energy 
index between p {2 — ?>), depending upon the relative positions of the breaks in its spectraP^. 
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The observed rate of relativistic SNe in the local universe is consistent with the required rate of 
X-Ra>^^ and radio (see Suppl. Info) transients accompanying the UHECR accelerators. 

It has been found that the arrival direction of the Auger events correlate well with the lo- 
cations of nearby AGN^, this suggests that they come from either AGNs or objects with similar 
spatial distribution as AGNs. Note that the HiRes eventPSIdo not show such a correlation. Further- 
more, UHECRs correlate welP^ with the locations of neutral hydrogen (HI) rich galaxies from the 
HI Parkes All Sky Survey (HIPASS). Our proposal relies on the acceleration of UHECRs in the 
mildly relativistic outflow from a subset of SNe Ibc, for which we have determined the size and 
magnetic field evolution using our radio observations, rather than hypothetical magnetars, with as 
yet unknown magnetic fields, supposedly formed during sub-energetic GRBs^^ considered in that 
work. SNe Ibc occur mostly in gas rich star forming spirals. In particular the 21 cm fluxes of 
NGC3278 (hosting SN 2009bb) obtained from the HyperLeda database amount to ~ 1.9 x lO^M© 
of HI. Hence, the observed correlation of UHECR arrival directions with HI selected galaxie^ is 
consistent with our hypothesis. 

In this letter we have shown that the newly established subset of nearby SNe Ibc, with engine- 
driven mildly relativistic outflows detected as sub-energetic GRBs, XRFs or solely via their strong 
radio emission, can be a source of UHECRs with energies beyond the GZK limit. Our study 
demonstrates for the first time, a new class of objects, which satisfy the constraints which any 
proposed accelerator of UHECRs has to satisfy. If SN 2009bb is characteristic of this newly 
discovered class, a radio survey to detect all such events should be undertaken (see Suppl. Info). 
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As an example, an all sky radio survey atu — 1 GHz with a sensitivity of 1 mJy and cadence of 2 
months, can detect all such transient sources which can accelerate Neon nuclei to 60 EeV, within 
200 Mpc of the Earth (~ 4 per yr). Such a survey will also detect radio emission from ordinary 
SNe and non-relativistic transients. However, their faster rise to peak will easily separate out the 
relativistic SNe for multi-frequency follow up. 
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Figure 1: Hillas Diagram: Mildly relativistic sources {(3/V ~ 1) must lie above the solid red line, 
to be able to accelerate Iron nuclei to 60 EeV by diffusive shock acceleratiorPl, according to Ez ^ 
f3eZBR/T^. In comparison, non-relativistic SNe {13 /T ~ 0.05) must lie above the dashed blue 
line to reach the same energies. Radius and magnetic field of SN 2009bb (red crosses, at 5 epochs, 
determined here from radio observations with VLA and GMRT assuming equipartition) and XRF 
0602 ISf^ (magenta ball) lie above the solid red line. Othei'^' radio SNe with SSA fits are plotted as 
blue balls. Only the SN 1993J magnetic fields are obtained without assuming equipartition'^. Note 
that all of the observed non-relativistic SNe (blue balls) including SN 1993 J (green balls) lie below 
the dashed blue line and are unable to produce UHECRs unlike the mildly relativistic SN 2009bb 
and XRF 060218 which lie above the red line. Sizes of crosses are twice the standard errors, sizes 
of balls are bigger than the standard errors. 
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Supplementary Information 



Radius-Magnetic Field Evolution: With a robust set of assumptions for the electron energy dis- 
tribution and magnetic field^, the radius of the forward shock wave at the time of the synchrotron 
self-absorption peak can be written ad^ 



i? ^ 4.0 X lO^^a-i/^M n I}^] ^] / ^ ^ ^ 



cm, (1) 



-5 GHz. 

where a = eg/e^ is the ratio of relativistic electron energy density to magnetic energy density, / 
is the fraction of the spherical volume occupied by the radio emitting region, Fop is the observed 
peak flux, and D is the distance. Using the same variables, the magnetic field is given by 

/ f \ -4/19 / p \ -2/19 / n \ -^/19 . „ . 

«---"'(i5) (^) {£} (5(k)°- 
The radio spectrum of SN 2009bb at all epochs from discovery paper (Fig. 2 of rei^ and this work, 
as obtained from observations using the Very Large Array (VLA) and the Giant Metrewave Radio 
Telescope (GMRT), is well fit by the SSA model, giving us a rare opportunity to explicitly measure 
the size and magnetic field of a candidate accelerator, instead of indirect arguments connecting 
luminosity with the Poynting flux. 

Equipartition: As already stated, non-relativistic SN have sizes and magnetic fields which are 
characteristically inadequate to accelerate charged particles to the highest energies. However the 
inferred magnetic fields are based mostly on equipartition arguments. The only young SN where 



the magnetic field was determined independent of the equipartition assumption was SN 1993J (this 
was however a type lib SN, unlike the type Ibc's being considered here). However, the magnetic 
field determined using a synchrotron cooling break was found to be ~ 9 times larger than the 
equipartition valueP'. This only helps by placing the SN 1993 J much closer to the Hillas line 
(but still below it) due to enhanced BR product. The energy requirement to explain the radio 
emission has a minimurri^S at the assumed equipartition factor of a = 1. In the absence of an 
independent measurement of the magnetic field, as in SN 1993^^, Very Long Base Interferometry 
of the outflow can constrain the deviation from equipartition using Equation ([T]). However, because 
of the very slow dependence of the radius on a, even for a slight deviation from equipartition in 
SN 2009bb, it can easily inflate the energy in the radio emitting plasma. In the case of SN 2009bb, 
even a more conservative assumption of equipartition, aided by its demonstrated mildly relativistic 
outflow enables it to be in a class of SN which can readily account for possible accelerators of 
UHECRs. 

Energy Budget in SN 2009bb: The radius evolution of SN 2009bb, as inferred from our radio 
observations (Table 1), is consistent with almost free expansion. This can only be explained if the 
mass of the relativistic ejecta is still much larger than the swept up mass. The computed Lorentz 
factor has barely decreased from 1.32 to 1.23 between days 20 and 222 post explosion. Using the 
model for coUisional slowdown of the ejecta, we modify Equation 1 15 of Rel'^^, to give 
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(71 - 1)^/^(71 + 1)^/' / (y - + ir'^'di, (3) 



m{Ri) + Mo j^^ 
where m(i?i) and m(i?2) are the swept up mass at the two respective epochs. We have neglected 

radiative losses, as they are unlikely to be important for protons in the time range of interest. 
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Moreover radiative losses would only help increase our initial energy budget. Performing the 
integral numerically from 71 to 72, the Lorentz factors at the two epochs and substituting for m{R) 
using the progenitor mass loss rate!^, we solve for the ejecta mass to get Mq ^ 1.4 x lO^^M©. 
Most of the mass in the relativistic outflow is due to baryons. The energy associated with these 
relativistic protons and nuclei is found to be Esaryons ^ 3.3 x 10^^ ergs. Compared to this blast- 
wave calorimetric value, the equipartition energy in the electrons and magnetic fields determined 
from SSA fit to the radio spectrum was was reportecf*^ to be E^q ^ 1.3 x 10"^^ ergs. This gives the 
electrons only a fraction e = ^ ^ 0.002 of the energy in the relativistic baryons. If EBaryons is 
distributed equally over 10 decades in energy, it can account for ~ 0.33 x 10^^ ergs of energy in 
UHECRs per logarithmic energy interval. Given the rate of relativistic SNe in the local universe, 
this is consistent with the volumetric energy injection rate for UHECRs. 

Rate of X-Ray transients: The number density of X-Ray flares associated with UHECR accelera- 
tors has been prescribe(i^ assuming that, the accelerated electrons have the same initial power-law 
index for their energy spectrum as the protons, and that they lose all their energy radiatively. Using 
the values of the physical parameters, motivated by SN2009bb,the number density of active X-ray 
flares with a luminosity > z/L^, is then given by recasting Equation 2.5 of Ref^to give 

fiAt ^ 3 X 10-^ { ( —r, \ ( ' Mpc-^ (4) 

V0.OO2/ V104^ergMpc"3yr~V VlO ergs-iy ^ 

SN 2009bb was observed with the Chandra ACIS-S instrument, at age 3 1 days, in the energy range 
0.3-10 keV. It had an X-ray luminosit>ii3 of Lx = '^A ± 0.9 x 10^*^ erg s"^ This luminosity 
and the rate of the relativistic SNe, considered in this work, together can account for the UHECR 
flux, if they remain active accelerators for At of order ~ 1 year. This is consistent with our radio 
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observations, which confinn that the BR product remains above the threshold throughout the ~ 1 
year of observation (See Table 1). 



Rate of radio transients: The required rate of radio transients, which can supply the observed 
Tinj is given by ri = Tinj/EsN- Assuming that the electrons and magnetic fields together have 
a fraction e of the energy of the relativistic protons (which is assumed to be divided equally into 
~ 10 logarithmic bins, assuming p ^ 2 for the protons), we compute the minimum required rate 
of such transients with peak radio luminosity Lop, which remain mildly relativistic at least until the 
SSA peak frequency drops to u, as 

/ r \ / ^ N / T \ -23/19 

—7 / 1 / \ I -'-'op 



n — 3 X 10' 

' 10^4 erg Mpc-3yr-i y V0.002y V 10^9 ergs/sec/Hz ^ 

Here r] = 6obs/deq is the ratio between the observed angular radius and the one obtained by as- 
suming equipartition between electrons and magnetic field^. This criterion works for an electron 
energy index between p = 2 (with the cooling break shifted below the SSA peak) to p = 3 (with 
the cooling break above the observed radio frequencies), so as to explain the observed spectral in- 
dex of ^ 1 in the optically thin part of the radio spectrum. Here, Equation (|5]) is the radio analogue 
of Equation 2.5 of Re^ (which is for X-ray transients). Hence, the observed rate of relativistic 
SNe can easily explain the energy injection rate. 

Survey Parameters: The prototypical mildly relativistic SN 2009bb has been discovered^ in a 
dedicated radio follow-up of type Ibc SNe. To firmly establish the rate of occurrence of such 
relativistic SNe in the local universe, a systematic large area radio survey is required. Here we 

iv 



estimate the maximum energy to which relativistic SNe can accelerate nuclei of charge Ze from 
Ez < ZeBR as /3/r ~ 1 for mildly relativistic outflows. Substituting the expressions for the 
radius (Equation [T]) and magnetic field (Equation [2l) we have 

/ f \ -5/19 / p \ 7/19 / n \ 14/19 

E^^QAxZa-'/'^-^] 1-^] I , , ) EeV, (6) 



0.5 J \mJj J V200 Mpc^ 
which is independent of the observed SSA peak frequency u. Further, assuming that R ~ T(3ct 
and we get the time to reach the synchrotron peak is 

Hence relativistic SNe will have faster rise times than non-relativist radio transients, allowing them 
to be easily identified for multi-frequency follow up with targeted observations. These considera- 
tions determine the cadence and sensitivity of the proposed radio survey as mentioned in the main 
text. 

Propagation and Survival of Nuclei: In the particle acceleration scheme outlined in this work, 
the highest energy particles are likely to be nuclei rather than protons. This is borne out by the 
latest Auger data which favors an increasing average rest mass of primary cosmic ray particles 
at the highest energies^*^. As for protons, the flux of ultra high energy nuclei are also suppressed 
over cosmological distances via interaction with background radiations. CMB photons appear 
has high energy 7-rays when Lorentz boosted into the rest frames of ultra high energy protons or 
nuclei. Protons above ~ 60 EeV, interact with CMB photons via the A resonance (7cmb + P ^ 
A+ — p + 7r° or 7cMB + P — )■ A+ n + 7r+) and give rise to the GZK limiP"^. Similarly ultra 
high energy nuclei can be photo-disintegrated by Lorentz boosted cosmic infrared background 



photons interacting mainly via Giant Dipole Resonances. The distance over which this effect 
suppresses the flux of ukra high energy nuclei is a function of the nuclear species and its energy. 
Detailed calculations^^ using updated photo-disintegration cross-sections indicate that the energy 
loss lengths for 100 EeV intermediate mass nuclei such as Ne, Si and Ca are around ~ 100 Mpc. 
As discussed in the main text, relativistic supernovae can provide enough number of UHECR 
sources within this distance, to be consistent with observation of independent arrival directions for 
the UHECRs. 

High Energy Particle Detection: UHECRs from the same source but with different energies will 
travel by different trajectories due to deflections by magnetic fieldP. For current estimates of 
the average intergalactic magnetic field, the mean delay in the arrival time of UHECRs, when 
compared to photons is found^ {Tdeiay) ~ 10^ yrs. Hence, barring chance coincidences, detected 
cosmic rays will not point back to known astrophysical transient^. However, detected UHECRs 
should point back (within the errors from deflection) to the host galaxies. As type Ibc supernovae 
occur mostly in HI rich spirals, the detected correlation with HI selected galaxie^ is consistent 
with our hypothesis. Similarly direct detection of UHECRs from say SN 2009bb is unlikely. 
However photo-hadron interaction between accelerated protons or nuclei and optical photons from 
the underlying SN may produce pions which then decay (7r+ — + Ue + + or — t- 
e" + i^e + ^'/x + ^v) to give high energy neutrinos. Neutrinos have no electric charge, hence they are 
not deflected by the intergalactic magnetic fields. They have very low rest masses compared to their 
very high energies and will travel at nearly the speed of light. These neutrinos will not be coincident 
with the initial core collapse as the number of accelerated charged particles which are the source of 
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neutrinos grows with time. The peak of the high energy neutrino flux will approximately coincide 
with the peak in bolometric luminosity (at around a week after explosion for SN 2009bb) as the 
most number of photons will be available for interaction with the accelerated hadrons. Hence, high 
energy neutrinos may by found in future by neutrino detectors like the IceCube, in directional and 
rough temporal coincidence with relativistic supernovae. 

Alternative Sources: AGNs The arrival directions of 20 of the 27 highest energy cosmic rays 
detected by The Pierre Auger Cosmic Ray Observatory, based in the southern hemisphere, were 
found to be within 3.2° of AGNs within 75 MpS^. This leads The Pierre AUGER Collaboration to 
conclude that they possibly come from either AGNs or objects with a similar spatial distribution. 
Yet no significant correlation is seen for the UHECRs detected by the HiRes stereo events and 
AGN^ in the northern hemisphere. Even, the degree of correlation in the Auger events now 
appears to be weakei® than what was seen by the earlier data. It has been suggested that these 
cosmic rays may be accelerated in the relativistic outflows from powerful AGNs. However, particle 
acceleration to such high energies {E = 10^° x E20 eV) in turbulent shocks with bulk Lorentz factor 
r would be accompanied by a minimum power lost to the Poynting flu?P® of 

L > lO^^r^E^o erg s-\ (8) 

Continuous sources of such luminosity would be easily detected within 200 Mpc and their absence 
rules out continuous AGN jets as the sources of a significant fraction of the UHECRs. Instead, a 
new class of very intense, short-duration AGN flares were proposed as possible sourceP. However, 
no such flare has been observed until now. 
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Alternative Sources: Classical GRBs Classical GRBs have also been suggested as one of the 



be accelerated by the Fermi mechanism in an ultra-relativistic outflow. For an astrophysical source 
driving a magnetized plasma outflow with a characteristic magnetic field B, at a velocity v = (5c 
(bulk Lorentz factor F), out to a radius R, the maximum energy to which a proton of charge e can 



GRBs satisfy the minimum luminosity criterion (Equation[8]) derived from the Poynting flux carried 
out by this outfloW^*^. However, most classical GRBs are found at cosmological distances with a 
mean redshift of 2.8 for those discovered by the SwifP^, hence most GRBs cannot contribute to the 
flux of cosmic rays above the GZK limit. If classical GRBs are indeed the source of the observed 
flux of the UHECRs, then the observed local rate of GRBs implies that each GRB is required to 
provide of orders of magnitude more energ)^', than what is available from a coUapsar scenaricP'. 

Alternative Sources: Hypernovae Hypemovae with a continuous ejecta profile (with Ek oc 
(r/3)~^) between the non-relativistic and relativistic material have been suggested as sources of 
UHECR^. In such a model each shell of different velocity accelerates particles upto a different 
energy and adds up to a final power law energy spectrum with slope of — 1 for E'^dN/ dE, which 
was claimed fits the observed UHECR spectra^. However, the observed spectrum of UHECRs is 
suppressed beyond the GZK limit''^ via interaction with the CMB photons and has to be corrected 
for the propagation effects to get the original injection spectrum. Hence Hypemovae with contin- 
uous ejecta profiles cannot reproduce flat injection spectrum of UHECRs which requires equal en- 
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most promising candidates for producing the highest energy cosmic ray^^, where protons would 



be accelerated by diffusive shock acceleratioii^ is given 




(9) 



ergies in each logarithmic energy bin. Galactic trans-relativistic SNe have also been considerecP^. 
However, this requires at least one trans-relativistic SNe per Galactic confinement time for the cos- 
mic ray energies being considered. For particles with energies beyond the GZK limit, the magnetic 
rigidities are so high, that their confinement time is comparable to the light crossing timeP of 10^ 
years^^. There are around ~ 100 SNe in this time, of which say upto 10 are SNe Ibc. Given the 
fraction of Ibc SNe which have relativistic outflow^, there are ~ 0.07 trans-relativistic SNe in 
this time. Clearly, if these events are galactic, the rate of such objects is too low and a galactic 
origin would neither explain the independent arrival direction^ nor the GZK suppression of the 
spectrcP^. 
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Observation 


Age 


op 




R 


B 


Ep 




Date (2009) 


(Days) 


(mJy) 


(GHz) 


(lO^^cm) 


(mG) 


(EeV) 


(EeV) 


05 April 


17 


>24.53 








>6.4 


>166 


08 April 


20 


17.87±0.95 


7.63±0.63 


34±3 


570±48 


5.7±0.1 


148±3 


10 IVlay 


52 


13.69±0.79 


3.33±0.17 


68±4 


256±14 


5.2±0.1 


134±3 


08 June 


81 


10.82±0.34 


1.93±0.07 


106±4 


152±5 


4.7±0.1 


123±1 


10 August 


144 


9.82±0.65 


0.90±0.06 


216±16 


72±5 


4.6±0.1 


119±3 


27 October 


222 


8.35±0.59 


0.53±0.04 


337±28 


43±3 


4.3±0.1 


112±3 



Table 1: Radius-IVIagnetic Field Evolution: Peak fluxes and peal< frequencies of SN 
2009bb are determined from VLA and GIVIRT observations by fitting a SSA spectrum 
to the observed fluxes. Fluxes until August can be found in Supplementary Info, of Re^^. 
The fluxes around 27 October 2009 are from new VLA (1.6 ± 0.1 mJy at 8.46 GHz and 
3.7 ± 0.2 mJy at 4.86 GHz) and GMRT observations (4.4 ± 0.3 mJy at 1 .28 GHz, 8.7 ± 0.7 
at 617 MHz and 5.8 ± 0.7 mJy at 332 MHz). Radius and magnetic fields are determined 
using Equations (fT|[2]> . Maximum energies to which protons and Iron nuclei can be accel- 
erated are computed using Equation dD, ± are standard errors. Note that both Ep and 
Epe decrease slowly by only ~ 24% in a span of ~ 200 days. 



